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ABSTRACT 

We present a new diagnostic diagram for mid-infrared spectra of infrared galaxies based on the equivalent 
width of the 6.2 /im PAH emission feature and the strength of the 9.7 /im silicate feature. Based on the position 
in this diagram we classify galaxies into 9 classes ranging from continuum-dominated AGN hot dust spectra and 
PAH-dominated starburst spectra to absorption-dominated spectra of deeply obscured galactic nuclei. We find 
that galaxies are systematically distributed along two distinct branches: one of AGN and starburst-dominated 
spectra and one of deeply obscured nuclei and starburst-dominated spectra. The separation into two branches 
likely reflects a fundamental difference in the dust geometry in the two sets of sources: clumpy versus non- 
clumpy obscuration. Spectra of ULIRGs are found along the full length of both branches, reflecting the diverse 
nature of the ULIRG family. 

Subject headings: galaxies: ISM — infrared: galaxies — galaxies: active — galaxies: starburst 



Lutz et al 



Dale et al 



1. INTRODUCTION 

Over the last decade several diagnostic diagrams have been 
proposed to quantify the contribution of star formation and 
AGN activity to the infrared luminosity of infrared galax- 
ies based on mid-infrared (to far-infrared) continuum slope, 
PAH line-to-continuum ratio, PAH to far-infrared luminos- 
ity ratio and the ratio of a hi gh to a low ionization forbid- 
den li n e such as [Nevl/[Ne ii l ?Genzel et alj|1998t 
1998; lLaurent et all l2000t iPeeters et alJ 120041: 
2006; Stur m et al.ll2006l) . However, none of these diagrams 
takes into account the effects of strong obscuration of the nu- 
clear power source. 

With the adve nt of the Infrared Spectrograph (IRS; 
iHouck et alj 120041) on board the Spitzer Space Telescope 
(IWerner et all 2004) astronomers have been handed a power- 
ful tool to study the 5-37 /im range for a wide range of galaxy 
types at an unprecedented sensitivity. This enables for the 
first time a systematic study of a large number of galaxies 
over the wavelength range in which amorphous silicate grains 
have strong opacity peaks due to the Si-O stretching and the 
O-Si-O bending modes centered at 9.7 and 18 /mi, respec- 
tively. 

Here we will introduce the strength of the 9.7 pm silicate 
feature as a tool to distinguish between different dust ge- 
ometries in the central regions of (ultra-)(luminous) infrared 
galaxies as part of a new diagnostic diagram and mid-infrared 
galaxy classification scheme. 

2. OBSERVATIONS AND DATA REDUCTION 

The results presented in this paper are based on both Spitzer 
IRS and ISO-SWS observations. The core sample is formed by 
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the IRS GTO ULIRG sample (PID 105; J.R. Houck PI), which 
comprises ^100 ULIRGs in the redshift range 0.02 < z < 
0.93. This sample is compared to samples of A GN and star- 
burst t emplates from the IRS GTO programs 14 (Brandl et al. 
I2006t IWeedmanetaT] 120051) and 96 and to selecte d AGNs 
from the GTO programs 82 and 86 dShi et al.ll2006h . Addi- 
tional IRGs, LIR Gs and ULIRGs we re taken from program 
159 (NGC 1377: TRoussel et al.ll2006h and from the Spitzer 
DDT program 1096. Archival ISO - SWS spectra of starburst 
nuclei were taken from Stur m et all (l2000h . 

The Spitzer observations were made with the Short-Low 
(SL) and Long-Low (LL) modules of the IRS. The spectra 
were extracted from DROOP level images provided by the 
Spitzer Science Center (34 using pipeline version SI 1.0.2, 
131 using pipeline version S14.0) and background-subtracted 
by differencing the first and second order apertures. The spec- 
tra were calibrated using the IRS standard stars HD 173511 
(5.2-19.5 pm) and £Dra (19.5-38.5 /im). After flux calibra- 
tion the orders were stitched to LL order 1, requiring order- 
to-order scaling adjustments of typically less than 10%. 

3. ANALYSIS 

For all the spectra in our sample we have measured the 
equivalent width of the 6.2 /im PAH emission feature as well 
as the strength of the 9.7 pm silicate feature and plotted the 
two quantities in a diagnostic diagram (Fig.|TJ>. 

The flux in the 6.2 /im PAH emission band is measured by 
integrating the flux above a spline interpolated local contin- 
uum from 5.95-6.55 pm. The equivalent width (EW) of the 
PAH feature is then obtained by dividing the integrated PAH 
flux by the interpolated continuum flux density below the peak 
(~ 6.22 /im) of the PAH feature. 

The apparent strength of the 9.7 /tm silicate feature is in- 
ferred by adopting a local mid-infrared continuum and eval- 
uating the ratio of observed flux density (f b s ) to continuum 
flux density (f cont ) at 9.7 /im and defining 

/ obs (9.7/xm) 
/cont(9.7/yw) 

For sources with a silicate absorption feature S s u can be inter- 
preted as the negative of the apparent silicate optical depth. 

Given the great diversity among our mid-infrared galaxy 
spectra, there is no "one-size-fits-all" procedure to define 
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FIG. 1 . — Diagnostic plot of the equivalent width of the 6.2 PAH emission feature versus the 9.7 fim silicate strength. Upper and lower limits are denoted 
by arrows. The galaxy spectra are classified into 9 classes (identified by 9 shaded rectangles) based on their position in this plot. Colors are used to distinguish 
the various classes. From class 1A to IB, 1C, 2A, 2B, 2C, 3A and 3B the colors used are: green, cyan, dark blue, yellow, pink, purple, red and orange. After 
class assignment the PAH equivalent widths were corrected for the effect of 6 fim water ice absorption on the 6.2 continuum. The extent of the individual 
corrections are indicated by dotted horizontal lines. The two dotted black lines are mixing lines between the spectrum of the deeply obscured nucleus of NGC4418 
and the starburst nuclei of M82 and NGC7714, respectively. Galaxy types are distinguished by their plotting symbol: filled circles: ULIRGs and HyLIRGs; 
filled triangles: starburst galaxies; filled squares: Seyferts and QSOs; filled diamonds: other infrared galaxies. 



the local continuum in all spectra. We therefore developed 
separate methods for three distinctly different types of mid- 
infrared galaxy spectra: continuum-dominated spectra, PAH- 
dominated spectra and absorption-dominated spectra. These 
methods are illustrated in Fig. [2] and described in its caption. 
We are forcing the local continuum to touch down at 14.0- 
14.5 /im, because the ISM dust cross-section decreases be- 
tween the two silicate peaks in this region. Detailed radiative 
transfer calculations verify that our interpolation procedure 
properly reproduces the emission that would be generated by 
dust stripped of its silicate features (Sirocky et al. in prepara- 
tion). 

3.1. Galaxy classification 

The galaxy spectra in Fig.Q] are classified into 9 different 
classes based on their 6.2 /im PAH EW and 9.7 /jm silicate 
strength. The parameter space covered by the various classes 
is indicated by shaded rectangles in Fig.Q] Average mid- 
infrared spectra for the eight populated classes are shown in 
Fig-12] The average spectra were constructed by normaliz- 
ing all spectra to unity at 14.5 /im flux before the averaging 
process. In order to maximize the signal-to-noise (S/N) of 
the average spectra, low S/N spectra were discarded from the 



process. Below we describe the 8 average spectra in the order 
(bottom to top) they are presented in Fig. [3] The ninth class, 
3C, is not populated. 

The class 1A spectrum is characterized by a nearly fea- 
tureless hot dust continuum with a very weak silicate absorp- 
tion feature at 9.7 /im. The class IB spectrum differs from 
the class 1 A spectrum by clearly showing the family of PAH 
emission features at 6.2, 7.7, 8.6, 11.2, 12.7 and 17.3 /im on 
top of a hot dust continuum. This hot dust continuum is nearly 
absent in the class 1C spectrum, allowing the PAH emission 
features to dominate the mid-infrared spectral appearance. 
Silicate absorption at 9.7 fim becomes noticable in the class 
2C spectrum as an increased depth of the depression between 
the 6-9 and 11-13/im PAH emission complexes. Another 
marked difference with the class 1C spectrum is the steep- 
ening of the 20-30 fim continuum and the appearance of an 
18 fim silicate absorption feature. In the class 2B spectrum the 
PAH features appear weaker than in the class 2C spectrum. In 
addition, the spectrum starts to show a 6 fim water ice and a 
6.85 fim aliphatic hydrocarbon absorption band. In the class 
3B spectrum these absorption features reach their maximum 
depths, while the 9.7 and 18 fim silicate features continue to 
increase in depth up to class 3A. Absorption features of crys- 
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FIG. 2. — Comparison of several methods used to determine the apparent 
9.7 fim silicate optical depth. The three methods differ by the way the local 
mid-infrared continuum is defined. For continuum-dominated spectra (top) 
simple logarithmic linear fits are made to the 5.0-7.5 fim and 26.1-31.5 (im 
Jy-spectrum. These are spline interpolated over the 7.5-26 fim range, using 
a continuum point at 14.0 fim to guide the spline in the inter-silicate-feature 
range. Pivots and fit are shown in green. For PAH-dominated spectra (2nd 
and 3rd spectrum) the 5.5-14.5 fim continuum is estimated from a power 
law interpolation over this range. From 14.5 to 26.0 fim this power law is 
continued as a spline function with a curvature determined by continuum 
pivots at 26.0 and 30.0 fim. Pivots and fit are shown in blue. For absorption- 
dominated spectra (4th and 5th spectrum) the local 5.6-26 fim continuum is 
determined from a spline interpolation to continuum pivots at 5.2, 5.6, 14.0 
and 26.0 fim. An additional continuum pivot at 7.8 fim is included for spectra 
with very little PAH emission (e.g. bottom spectrum). Pivots and fit are 
shown in red. In all three cases the 9.7 fim silicate strength is computed by 
evaluating the ratio of observed and interpolated continuum at its peak. From 
top to bottom, the five spectra shown are the IRS low-resolution spectra of 
PG 1351+640, NGC7714, NGC3628, UGC510I and NGC44I8. 



talline silicates (Sp oon et al.l20 06) appear in the spectra of the 
classes 3B and 3 A at 16, 19 and 23 /mi. Equivalent widths of 
PAH emission features and emission lines decrease from class 
2B to 3B and 3A. Note especially the change in shape of the 
7.7 fj,m PAH feature as it first broadens and then disappears 
going from class 2B to 3B and 3A. Finally, the class 2A spec- 
trum differs from the class 3A spectrum mainly by a clearly 
lower apparent depth of the 9.7 and 18 /mi silicate features. 
Note that the spectral structure in individual spectra may dif- 
fer substantially from the average properties of the classes. In 
Fig.[3]this is represented by the 1-er dispersion ranges around 
the individual average spectra. 

3.2. Galaxy distribution 

The galaxies shown in Fig. Q] are color-coded according to 
their galaxy classification. However, their positions may dif- 
fer from their original assignments, as in Fig.Q] the 6.2 /im 
PAH EW has been corrected for the effect of 6 /im water ice 



absorption on the 6.2 /im continuum through substitution of 
the observed 6.2 /im continuum by the 6.2 /im continuum that 
was defined to infer the 9.7 /im silicate strength. The ice- 
correction is justified assuming that the ice is part of the ob- 
scuring medium behind the PAH emitting region. The cor- 
rection is applied only to spectra clearly showing the imprint 
of water ice absorption, as found among the spectra of classes 
2B, 2C, 3A and 3B. The extent of the correction for individual 
spectra is indicated in Fig.Q]by horizontal dotted lines. 

The ULIRGs, Seyferts, quasars and starburst galaxies 
shown in Fig.Q] are not distributed randomly through the 
diagram. Instead, the galaxies appear to be distributed 
along two branches. One extending horizontally from the 
positions of the pro totypical Seyfert-1 nucleus NGC4151 
(Weed man et alj|2005l) to the prototypical starburst nucleus 
NGC7714 (iBrandlet all 120041) . the other branch extending 
diagonally from the prototypi cal deeply obscured nucleus of 
NGC 44 1 8 dSpoon et al.ll200 lb to the starburst nuclei of M 82 
dSturm et all 120001) and NGC7714. Very few sources are 
found above and to the right of the diagonal branch and in 
between the two branches (the class 2A domain). 

The three extremes distinguish between AGN-heated hot 
dust spectra (class 1A), PAH-dominated spectra (classes 1C 
and 2C) and absorption-dominated spectra (class 3A). The 
galaxies lined up along the two branches in between the 
extremes show signatures of both extremes to varying pro- 
portions. The galaxies along the horizontal branch may be 
thought of as combinations of AGN and starburst activity, 
while the galaxies along the diagonal branch may be thought 
of as intermediate stages in between a fully obscured galac- 
tic nucleus and an unobscured nuclear starburst. The latter is 
illustrated by the two dotted lines in Fig.Q] which represent 
mixing lines between the spectrum of NGC 44 18 on one ex- 
treme and either NGC 7714 or M 82 on the other. The large 
majority of galaxies on the diagonal branch fall in between 
these mixing lines. One notable exception is Arp220. 

Relatively few galaxies are found in the section between 
the two branches: the domain of the class 2A sources. Spec- 
troscopically, there are several simple scenarios which will 
populate this class from adjacent classes. A class 3A galaxy 
can be turned into a class 2A galaxy by adding a featureless 
hot dust continuum to the class 3 A absorbed continuum spec- 
trum. The main effect of this continuum dilution is the fill- 
ing in of the deep 9.7 /im silicate feature, resulting in a less 
pronounced feature. Increasing the dilution not only further 
decreases the depth of the silicate feature but also increases 
the continuum at 6.2 /im, resulting in both a less negative sil- 
icate strength and a lower 6.2 /im PAH EW. The effect of this 
type of continuum dilution is illustrated in Fig.Q]by a curved 
arrow. Alternatively, a class 1A galaxy may be turned into 
class 2A galaxy by passing its mid-infrared spectrum through 
a foreground cold dust screen. The path that a galaxy will 
take is straight up, since the PAH EW is not affected by cold 
foreground extinction. 

The presence of foreground dust or dust mixed in within 
the PAH emitting region may also explain the difference be- 
tween class 1C and 2C starburst galaxies, and the difference 
between the classification of the prototypical starburst galax- 
ies NGC 7714 and M82 in particular. In the latter case this 
may be primarily an orientation effect, as the nuclear starburst 
in NGC 7714 is observed face-on, whereas the starburst ring 
in M 82 is seen under an inclination of ^80 degrees. 

Finally, the complete absence of galaxy spectra classified as 
3C must imply that strong, pure foreground, cold dust obscu- 
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FIG. 3. — Average spectra for eight of the classes of galaxy spectra defined 
in Fig.ff] along with the 1-cr dispersion range (grey shading). Most of the 
dispersion is due to the overall slope of the spectrum. Spectra at the bottom 
of the figure (classes 1A and IB) are dominated by hot dust emission, spectra 
in the middle (classes 1C and 2C) by PAH emission and spectra towards 
the top by absorption features (classes 3B, 3A and 2A). Vertical dotted lines 
denote the central wavelengths of the family of PAH emission features at 6.2, 
7.7, 8.6, 11.2 and 12.7 /jm. 

ration is not common in starburst nuclei. The obscuring dust 
instead contributes to the continuum underlying the 6.2 fim 
PAH feature, reducing its EW. 

4. DISCUSSION AND CONCLUSIONS 

We have constructed a diagnostic diagram of two mid- 
infrared spectral quantities, the equivalent width of the 6.2 fim 
PAH feature and the apparent 9.7 fim silicate strength, for 
the purpose of classifying infrared galaxies according to their 
mid-infrared spectral shape (Fig.[TJ. 

The large majority of galaxies (>90%) are dispersed around 
two branches: a horizontal one extending from continuum- 
dominated spectra to PAH-dominated spectra and a diago- 
nal one spanned between absorption-dominated spectra and 
PAH-dominated spectra. Seyfert galaxies and quasars (square 
plotting symbols in Fig.[T|i are found exclusively on the hor- 
izontal branch characterized by the absence of pronounced 
silicate features. Starburst galaxy spectra (triangular plotting 



symbols) are concentrated toward the extreme right tip of the 
two branches, the locus of PAH-dominated spectra. Of the 
ULIRGs and HyLIRGs in our sample (round plotting sym- 
bols) about a quarter occupy the same space as the AGNs and 
starbursts in our sample. The large majority, however, are 
distributed along the diagonal branch, characterized by in- 
creasingly apparent silicate absorption and less pronounced 
PAH emission features, with deeply obscu red galactic nu- 
clei such as NGC4 418 (ISpoon et alJl200l and NGC1377 
dRoussel et al.l2006l) as the only non-ULIRG, non-Seyfert end 
members on this branch. Note that low-metallicity galaxies 
(not included in this study) would be distributed along the 
horizontal branch. 

The existence of two distinct regimes of apparent silicate 
absorption, as suggested by the two distinct branches in Fig.Q] 
is intriguing and likely points to clear differences in the nu- 
clear dust distribution between galaxies on the two branches. 
Clump y distributions produce only sh allow absorption fea- 
tures dNenkova. Ivezic & Elitzur 120 02). and are the likely ex- 
planation for the horizontal branch. The deep silicate absorp- 
tion found on the slanted branch requires the nuclear source 
to be deeply embedded in a smooth dist ribution of dust that i s 
both geometrically and optically thick dLevenson et al.f ;2006). 
The transition from clumpy geometry to one dominated by 
a smooth dust distribution is a possible explanation for the 
sources located between the branches. 

The distribution of ULIRG spectra along the full length of 
both the horizontal and the diagonal branch once again il- 
lustrates the diverse nature of the ULIRG family. It, how- 
ever, also raises the issue of ULIRG evolution. For instance, 
what evolutionary path has taken IRAS 08572+3915 to its ex- 
treme class 3A position in which the central power source is 
deeply buried? Did the interaction start with both nuclei clas- 
sified as class 1C/2C starburst galaxies, which then gradually 
moved up the diagonal branch as the interaction stengthened 
and more dust accumulated on the remnant nuclei? And how 
will IRAS 08572+3915 evolve from there, once the obscur- 
ing screen breaks up and the hidden power sources are re- 
vealed? In case of a hidden starburst the source likely will 
move back diagonally toward the starburst locus. However, if 
the dominant power source is an AGN, will the source cross 
into the sparsely populated class 2A regime, evolving directly 
toward the AGN-dominated class 1A regime, or will it first 
undergo a starburst before settling somewhere along the hor- 
izontal branch? The sparse population of the class 2A do- 
main either indicates that the crossing time is brief - with few 
ULIRGs caught in transformation - or that deeply obscured 
ULIRG nuclei mostly evolve to starbursts (first). 

In a more extensive paper (Spoon et al. in preparation) we 
will present correlations of additional parameters within the 
diagram (e.g. mid-IR line ratios, Lir, Lpah/Lir, absorption 
features, nuclear separation, optical classification). 
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